Condensate fraction in a 2D Bose gas measured across the Mott-insulator transition 
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We realize a single-band 2D Bose-Hubbard system with Rb atoms in an optical lattice and measure 
the condensate fraction as a function of lattice depth, crossing from the superfluid to the Mott- 
insulating phase. We quantitatively identify the location of the superfluid to normal transition by 
observing when the condensed fraction vanishes. Our measurement agrees with recent quantum 
Monte Carlo calculations for a finite-sized 2D system to within experimental uncertainty. 



Measurements of condensed matter systems realized 
by cold atoms in optical lattices are now performed 
with sufficient accuracy to compare with ab-initio cal- 
culations 0, 0, 0|. Bosonic atoms in a lattice nearly 
perfectly realize the iconic Bose-Hubbard (BH) Hamil- 
tonian. Here, we study the system's momentum distri- 
bution, measure the condensate fraction, and accurately 
identify the transition point from the the low tempera- 
ture superfluid (SF) phase by identifying when the con- 
densate fraction vanishes. 

The SF to Mott insulator (MI) transition can be ac- 
cessed by changing the depth of the optical potential Q , 
and has been observed in ID @, 2D and 3D @. A 
range of studies have verified a detailed understanding 
of the MI phase in 2D and 3D [J [2, 0, Q- In contrast, 
the SF phase and the details of the transition to MI have 
gone largely unstudied. Indeed, the only quantitative 
measurement locating the transition is in 3D and is not 
in agreement with calculations Here we focus specif- 
ically on the superfluid phase of a 2D system and its 
transition to a normal state: we observe the expected 
increasing momentum spread and vanishing condensate 
fraction as the system leaves the SF phase. Our mea- 
sured transition point agrees with the best available cal- 
culations 0], thereby locating a point on the non-zero 
temperature 2D BH phase diagram. Interestingly, the 
condensate fraction in our non-zero temperature system 
vanishes more sharply than expected for a zero temper- 
ature inhomogenous system, confirming that the super- 
fluid regions are rap idl y driven normal as soon as an in- 
sulator appears 21, [22 1 . 

The physics of interacting systems frequently depends 
spectacularly on dimensionality: in 3D the SF is a 
conventional Bose-Einstein condensate (BEC); in 2D, a 
Berezinskii-Kosterlitz-Thouless (BKT) SF; finally, in ID 
there is no true SF. In contrast, only the detailed prop- 
erties of the MI phase depend on dimensionality. In the 
T > 0, 2D case studied here, the very existence of Bose 
condensation is a consequence of the finite size of our 
trapped system. We associate the presence of a bimodal 
momentum distribution with the SF phase, and use fits 
to the distribution to identify the Bosc-condcnscd frac- 
tion and thereby measure the transition point between 
SF and normal. 



At low temperature the transition from SF is to a nor- 
mal state which crosses over to a MI phase as the lattice 
depth increases 
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|22|. As a result any T > mea- 
surement based on condensate fraction will identify the 
SF to normal transition but be largely insensitive to the 
subsequent crossover into the MI phase. 

We study samples of ultra-cold rubidium atoms in a 
combined sinusoidal plus harmonic potential. For atom 
occupancy per lattice site larger than unity 0, [Io| , the 
low temperature SF phase (shallow lattice) is expected to 
evolve into a structure composed of alternating shells of 
SF and integer-occupied MI (deep lattice). As the lattice 
deepens, each successive MI region appears and grows, 
as probed in Ref. § . At T = the amount of SF varies 
smoothly with lattice depth giving no abrupt changes in 
the momentum distribution to indicate a phase transi- 
tion. In this work, we simplify the situation by work- 
ing near unit filling, where the only insulating phase is 
unit occupied MI; thus, any observed signature can only 
be the transition from SF to normal. Absent the lat- 
tice, recent experiments have shown that weak contact 
interactions lead to a decrease in the 2D BEC transition 
temperature [l3| • Lattice potentials increase the relative 
importance of interactions; indeed, the onset of the MI 
phase corresponds to driving the critical temperature to 
zero. 

The BH model describes lattice-bosons that have a 
hopping matrix clement t, and an on-site interaction en- 
ergy U. The physics of the BH model depends only on 
U/t [lj|. In an infinite, homogenous T = 2D system, 
the transition from SF to MI occurs at the critical value 
(U/t) c w 16.5 d Q3 [U El. Remarkably, we observe 
a sharp transition at U/t = 15.8(20) [H] in our T > 0, 
finite-sized, harmonically trapped system. 

Our data consists of images of atom density after sud- 
den release and time-of- flight (TOF), approximating the 
in-situ momentum distribution. Figure Q] shows 2D mo- 
mentum distribution (right), and cross-sections through 
each distribution (left). As evidenced by Fig. [TJa and 
-b, each diffraction order in the momentum distributions 
consists of a narrow peak on a broad pedestal. Fit- 
ting to a bimodal distribution (see below), we deter- 
mine /, the fractional contribution of the narrow com- 
ponent, and identify / as the "condensate" fraction. We 
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FIG. 1: Momentum distributions and cross-sections at U/t — 
4(1), 8(1), and 20(2). Each row shows a single momentum 
distribution normalized by the total atom number; the lines 
in the top-right panel indicate trajectories along which four 
cross-sections were taken. The left panel shows the average 
of these four sections (black solid line); the red-dashed lines 
denote the fit to the bimodal distribution. 



associate images with non-negligible / as being in the 
SF phase We emphasize, however, that superfluid- 
ity is a transport phenomena and cannot unambiguously 
be associated with features in the momentum distribu- 
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|22j. This association is also imperfect at 
T > because in our 2D trapped system we expect to 
observe a discernible "condensate" fraction even after the 
vortex-pairs of a BKT SF unbind [TlJ , destroying the 2D 
SF. / vanishes only when the resulting phase-fluctuating 
quasi-condensate vanishes [l2l [l3| . 

To characterize the transition from SF to normal, we 
extract two independent quantities from an analysis of 
TOF images: /, and an "energy scale" a. We also mea- 
sure a related quantity, the full width at half maximum 
(FWHM) r of the quasi-momentum distribution, which 
we compare to theory. As the lattice depth is increased 
we find that / vanishes concurrently with a sudden in- 
crease in r, abrupt signatures that we associate with the 
transition. 

We produce nearly pure 3D 87 Rb BECs with Nt = 



1.2(4) x 10 5 atoms in the |F = l,m F = -1) state 0- A 
pair of linearly polarized, A = 820 nm laser beams forms 
a 30(2)E R deep vertical optical lattice along z that di- 
vides the 3D BEC into about 70 2D systems (turn-on 
time = 200 ms). The single photon recoil wave- vector 
and energy are k R — 2tt/\ and Eji = h 2 k^ 2 /2m = 
h x 3.4 kHz; m is the atomic mass and h is Planck's 
constant. The largest 2D system, containing w 3000 
atoms, has a chemical potential /^2D = h x 600(100) Hz 
and we measure a temperature fc^T = ks x 33(4) nK = 
h x 700(70) Hz. Since the first vibrational spacing 
h x 33(1) kHz ^> /i2D ; IcbT, this system is well into the 
2D regime. In addition, a weaker, square 2D lattice in 
the x-y plane is produced by a second beam arranged in 
a folded-retroreflected configuration (23[, linearly polar- 
ized in the x-y plane (turn-on-time = 100 ms |30(). The 
intensities of both lattices follow exponentially increas- 
ing ramps, with 50 ms and 25 ms time constants respec- 
tively, and reach their peak values concurrently. These 
time-scales are chosen to be adiabatic with respect to 
mean-field interactions, vibrational excitations, and tun- 
neling within each 2D system. The final depth of the 
x-y lattice determines U jt and ranges from V = to 
25(2)Efi [3l|]. The lattice depths are calibrated by puls- 
ing the lattice for 3 fj,s and observing the resulting atom 
diffraction [24j ]. 

We calculate U / t using a 2D band-stucture model and 
the s-wave scattering length [14|. The uncertainty in U/t 
stems from the uncertainty in lattice depth [32|. The 
resulting uncertainty in U/t is ±10%. 

Once both lattices are at their final intensity, the 
atomic system consists of an array of 2D gasses each in a 
square lattice of depth V and with a typical density of 1 
atom per lattice site. The atoms are held for 30 ms, and 
all confining potentials are abruptly removed (the lattice 
and magnetic potentials turn off in < 1 p,s and ~ 300 /is, 
respectively). As a result, the initially confined states 
are projected onto free particle states which expand for 
a 20.1 ms TOF [33], when they are detected by resonant 
absorption imaging. Apart from effects of atomic interac- 
tions during expansion and the initial size of the sample, 
initial momentum maps into final position, so each im- 
age approximates the x-y projection of the momentum 
distribution. We fit the each momentum distribution to 
a simple function which describes the distributions over 
the full range of U/t studied here, with just three free 
parameters. 

First, we model the broad background as a ther- 
mal distribution of non-interacting classical particles 
in a single 2D sinusoidal band where states are la- 
beled by their quasi-momentum q x and q y , n(q x ,q y ) oc 
exp [2(cosirq x /kii + cos nq x /kn)/ a]; this contributes two 
fitting parameters: a and the non-condensed atom- 
number. In the shallow lattice limit, a gives the tem- 
perature, a = ksT/t. This fit does not distinguish 
atoms thermally occupying higher momentum states 
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from atoms occupying these states in the ground state 
wavefunction, i.e., from the quantum depletion of the 
SF. In fact, n(q xi q y ) multiplied by a suitable Wannicr 
function, correctly describes the momentum distribution 
of atoms in the MI phase to first order in t/U where a is 
unconnected to temperature, and is given by a — U/At. 
Our function fits the random phase approximation (RPA) 
momentum distribution fairly well even as higher order 
terms become important 0, 0] • 

The second portion of the momentum distribution 
consists of a narrow peak, which we interpret as Bose- 
condensed atoms. We therefore take the narrow peak to 
be the inverted parabola of a Thomas-Fermi profile (of 
fixed width for all comparable data 34]), characterized 
by a single fitting parameter, condensed number. 

The observed condensate peak width after TOF stems 
largely from initial system size, not interaction effects 
during TOF or the initial momentum spread. Here inter- 
actions during TOF are reduced due to rapid expansion 
along z after release from the tightly confining vertical 
lattice. Our analysis further reduces these interaction ef- 
fects by excluding data inside the 1 st Brillouin zone, with 
the highest density. This decreases the measured FWHM 
of the peak from 30(1) /<m to 22(1) /jm and the inferred 
momentum width from 0.26k R to 0.2lk R . Changing the 
TOF from 20.1 ms to 29.1 ms only increased the FWHM 
from 22(1) /im to 28(1) (im (decreasing the observed mo- 
mentum width from 0.21k R to 0.l7k R ). 

Figure [2]-a shows that as V increases, / vanishes at a 
critical value V cl -n, while the total atom number remains 
constant. We verified that this disappearance does not 
result from excessive irreversible heating of the system by 
exceeding V CI it, then lowering the lattice and observing a 
condensed fraction Q. 

To gain a qualitative understanding of the vanishing 
condensate fraction, we performed a non-zero tempera- 
ture mean-field theory (MFT) simulation of an array of 
2D BH systems in a 3D harmonic trap [13]. To model 
the non-zero temperature experimental system, we de- 
termine the entropy at small U/t that gives the observed 
« 45% condensate fraction, and assume this entropy is 
unchanged as V increases. The red-dashed line in Fig. [2} 
a shows the MFT condensate fraction vs. V at constant 
entropy. Given that T — MFT overestimates the tran- 
sition ((£//i)MFT = 23.3, compared to (U/t) c = 16.5 from 
more accurate calculations), the curve unexpectedly lies 
on the data. MFT also gives / as function of U/t in units 
of (U/t) c . We identify the transition point by fitting this 
function to the data allowing (U /t) c to vary, yielding 
(U/t) c = 15.8(20) (a lattice depth V ciit = 9.0(5) E R ). 

Figure [2]-b displays a from the uncondcnsed back- 
ground portion of the distribution. At large V we re- 
cover the behavior expected in the MI phase; this mea- 
surement is equivalent to observations of the modulated 
momentum distribution in the MI phase [H, 0| ■ At higher 
total atom number, our system would develop doubly 
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FIG. 2: Condensed fraction / and a vs. V (bottom axis) 
or U /t (top axis) . The dots denote values determined from 
2D fits to the full momentum distribution: small dots result 
from one image and the iarge dots indicate data averaged over 
about 20 separate images. The uncertainties are their RMS 
variation, and are indicative of the singie-image uncertain- 
ties, (a) Condensate fraction. The red-dashed fine is com- 
puted from our MFT model, (b) fit parameter a. There 
are two distinct regimes: at low U /t it is nearly constant 
(blue dashed line), from which we infer an initiai tempera- 
ture ksT « 2t; and at large U/t it monotonically increases, 
consistent with predictions of perturbation theory in the MI 
phase (red-dashed line). 



and triply occupied MI shells, expected to manifest as 
kinks in this curve, a is monotonic with V, varying 
smoothly across V cr it- This is in agreement with RPA 
theory where the onset of superfluidity affects only states 
near zero quasi-momentum. Figure ®b shows that when 
V < AE R (U/t < 3), k B T/t w 2.0(3). Extrapolating 
to V = gives k B T = E r o/-k 2 w k B x 33 nK (valid 
when T <C E R ). This temperature is well below the 
^ x 45 nK expected for non-interacting particles in our 
2D harmonic trap with / = 0.45, this reduction is similar 
to that observed in Ref. [13| . which focused on the crit- 
ical temperature in interacting 2D atomic systems with 
no 2D lattice. 

A related characterization of the system is the FWHM 
r of the full quasi-momentum distribution 0, 0, [l^, 
|25[. Figure [3] shows the width of the 2D distributions (see 
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FIG. 3: Quasi-momentum width vs. V (bottom axis) or U/t 
(top axis). The symbols denote the average FWHM of the 
quasi-momentum distribution along the axes of highest sym- 
metry (Top: averaged along x+y and x-y; Bottom: averaged 
along x and y). The small and large dots and uncertainties 
are as explained at Fig. The red-dashed line is the hori- 
zontally displaced RPA momentum width, as discussed in the 
text, and the vertical grey line denotes the location of the 
SF-normal transition identified from the sudden increase in 

r. 



Ref. 0) as a function of V . In the SF phase T hardly 
depends on V since the dominant feature of the distri- 
bution is the condensate peak. T only begins to change 
very close to the SF to normal phase transition when 
the heights of the two components of the bimodal dis- 
tribution become comparable when the condensate dis- 
appears, consistent with calculations in homogenous and 
trapped systems [9| . We calculate T in the MI phase from 
the RPA [2^ | quasi-momentum distribution which accu- 
rately describes the large U/t limit (y/2kR along x + y 
and kfi along x). In the RPA, V is a function of U/t in 
units of (U/t) c . The red dashed lines are fits to the mea- 
sured widths using two free parameters (joint between 
both panels): in the MI region we use the RPA functional 
form with (U /t) c as the first fit parameter, and the con- 
stant width in the SF phase is the second. We obtain 
(U/t) c = 16.7(20), in accord with the (U/t) c = 15.8(20) 
from our fit to the condensate fraction. 

We identify the point when the condensate fraction 
vanishes (and T abruptly increases) with the onset of 
the SF to normal transition, i.e., when a normal region 
begins to rapidly expand in our inhomogeneous system. 
(Our measured visibility, computed as in Refs. P, Q, 



abruptly drops from near unity at U/t 16.) In- 
creasingly accurate numerical calculations give values of 
(U/t) c : 16.25(10) and 16.77 [18]. Perhaps most rel- 
evant are QMC calculations which include the effects of 
harmonic confinement; in this case, Wessel et al. find 
that a MI region first forms at (U/t) c = 17.2 (the exact 
value of (U/t) c depends on the details of the harmonic 
potential). Both values lie within our experimental un- 
certainty. 

The calculations [3, [l8| are at zero temperature, and 
while they agree with our observed (U/t) c , they do not 
predict a sudden increase in peak width or a vanishing 
condensate fraction at (U/t) c . At T — and as U/t 
increases past (U/t) c , where an inhomogeneous system 
first develops a unit-occupied Mott core, the shell of SF 
persists to large U/t. Thus, at T = 0, / drops rapidly 
at (U/t) c , but does not vanish. Our system, however, is 
at small but non-zero temperature, with a reduced con- 
densate fraction of ~ 45% for small V. Our MFT model 
shows that this temperatures quickly drives the SF shells 
to the normal phase as U/t increases past (U/t) c . As 
a result the SF shells rapidly disappear as the normal, 
and then Mott regions form (see Refs. [HI S3, HI] ) . Tnat 
this feature is seen in preliminary non-zero temperature 
QMC calculations [IH underscores the need for further 
non-zero temperature calculations to compare with ex- 
periment. 

This experiment constitutes the measurement of a sin- 
gle point of the non-zero temperature 2D BH phase dia- 
gram. We expect future experiments will expand on this 
result at different temperatures, densities, in different di- 
mensions, and in traps with more homogenous density 
distributions; new theory should aide in the interpreta- 
tion of these experiments. 
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